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Background: The relativistic Hartree-Fock-Bogoliubov (RHFB) theory has recently been developed and it provides a unified 
and highly predictive description of both nuclear mean field and pairing correlations. Ground state properties of finite nuclei 
can accurately be reproduced without neglecting exchange (Fock) contributions. 

Purpose: Finite-temperature RHFB (FT-RHFB) theory has not yet been developed, leaving yet unknown its predictions for 
phase transitions and thermal excitations in both stable and weakly bound nuclei. 

Method: FT-RHFB equations are solved in a Dirac Woods-Saxon (DWS) basis considering two kinds of pairing interactions: 
finite or zero range. Such a model is appropriate for describing stable as well as loosely bound nuclei since the basis states have 
correct asymptotic behaviour for large spatial distributions. 

Results: Systematic FT-RH(F)B calculations are performed for several semi-magic isotopic/isotonic chains comparing the 
predictions of a large number of Lagrangians, among which are PKAl, PKOl and DD-ME2. It is found that the critical 
temperature for a pairing transition generally follows the rule = 0.60A(0) for a finite-range pairing force and Tc = 0.57A(0) 
for a contact pairing force, where A(0) is the pairing gap at zero temperature. Two types of pairing persistence are analysed: 
type I pairing persistence occurs in closed subshell nuclei while type II pairing persistence can occur in loosely bound nuclei 
strongly coupled to the continuum states. 

Conclusions: This first FT-RHFB calculation shows very interesting features of the pairing correlations at finite temperature 
and in finite systems such as pairing re-entrance and pairing persistence. 

PACS numbers: 05.70.Ce, 21.60.Jz, 24.10.Pa, 64.60.Ht 


I. INTRODUCTION 

The thermodynamical properties of excited nuclei have 
drawn over the last decades a renewed attention due to the 
advanced accurate measurements of level densities at low ex¬ 
citation energies M. Pairing correlations play an essential 
role in many Fermion systems and have thus a strong influ¬ 
ence on nuclear structure at low excitation energies Ma- HOr . 
Pairing correlations in finite systems such as nuclei or Wigner- 
Seitz cells, and in infinite ones such as in neutron star mat¬ 
ter, may exhibit different behaviors reflected in the specific 
heat and the level density lEIHIl. Moreover, the phase tran¬ 
sition is a complex and rich phenomenon, where pairing re¬ 
entrance in asymmetric matter, in odd-nuclei, rotating nuclei, 
and even in doubly magic nuclei close to the drip line may 
occur |l 3 [^[^[^. The interplay between temperature and 
shell effects in superfluid systems, giving rise to re-entrance 
or its opposite phenomenon — suppression, still remains to 
be studied. 

The competition between temperature and pairing correla¬ 
tions in nuclei at low excitation energies has been studied for 
several decades, with the pioneerii^ works based on finite- 
temperature BCS (FT-BCS) theory ifisll and finite-temperature 
Hartree-Fock-Bogoliubov (FT-HFB) theory ifl^ . It was pre¬ 
dicted that the critical temperature Tc for pair correlation 
quenching could be expressed, as in uniform matter, as a 
function of the average pairing gap at zero temperature A(0) 
following the relations: Tc - 0.57A(0) for a constant pairing 
force lllSn . The more evolved Bogoliubov approach has later 
been applied to finite nuclei confirming the existence of such 
relations between Tc and A(0). For the simplified degenerate 
model, the relation was found to be Tc - 0.50A(0) ifl^ . while 


in rare earth transition nuclei the interplay with deformation 
induces shape transitions in the superfluid phase, leading to 
the ratio 7’e/A(0) = 0.57 for protons and 0.63 for neutrons lEoll . 
In addition, pairing correlations are expected to play an im¬ 
portant role in the decay of compound nuclei formed in heavy 
ion collisions, as illustrated in the seminal work presented in 
Ref. im. More recently, the BCS and HFB approaches have 
been extended to self-consistent mean held models in order 
to improve the description of the pairing transition in spheri¬ 
cal nuclei ll23 - l^ . as well as in deformed nuclei where shape 
transitions have been predicted ll25l42^ . In summary, the ratio 
Tc/A{Q) lies in the interval 0.50 - 0.60, where the uncertainty 
originates mainly from the detailed level structure of spherical 
and deformed nuclei which depends itself on models. 

It is worth noticing that in most of the quoted studies, the 
calculations were performed either on the harmonic oscilla¬ 
tor basis, or within the non-relativistic framework. Due to 
the limitation of the harmonic oscillator basis in giving an 
appropriate asymptotic behavior of the single particle (s.p.) 
wave functions, the nucleon densities at large distance con¬ 
verge very slowly with respect to the size of the basis. The 
situation becomes even more serious in the weakly bound nu¬ 
clei close to the drip line SEl. Nowadays a realistic frame¬ 
work is to perform the calculations in an appropriate basis that 
can provide a reasonable description of both the overall and 
asymptotic behaviors of the density profiles, for instance the 
Woods-Saxon (WS) basis j^l^ . In some applications, the 
small component of the Dirac spinors were usually neglected 
in determining the relativistic Hartree-Bogoliubov (RHB) and 
finite-temperature RHB (FT-RHB) pairing tensor (| 3113 ] • We 

therefore present, in this paper, the first fully RHFB calcu¬ 
lations at finite temperature (FT-RHFB): both the large and 
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small components of the Dirac spinors contribute to the pair¬ 
ing channel, and the contributions of the Fock terms are natu¬ 
rally included. In addition, the Dirac WS (DWS) basis ll^ is 
employed to better describe weakly bound nuclei. 

This paper is organized as follows. The general formalism 
of the FT-RHFB theory is outlined in Sec. |II] In Sec. |III] we 
compare the results obtained with different covariant density 
functional (CDF) models |35[ -Miy using the finite- and zero- 
range pairing interactions 11^|4^ . Systematic FT-RHFB cal¬ 
culations are performed for several isotopes/isotones. More¬ 
over, the questions of pairing persistence and re-entrance phe¬ 
nomena with increasing temperature are addressed. Finally, 
the main conclusions are drawn in Sec. |IV] 


II. GENERAL EORMALISM AND NUMERICAL DETAILS 

We briefly recall here the general features of RHFB the¬ 
ory and thermodynamics, then their generalization to the finite 
temperature case. 


A. RHEB framework 

The relativistic Hartree-Fock (RHF) theory ll^ |4^ l43l] is 
designed to describe bulk and s.p. nuclear properties. The rel¬ 
evant degrees of freedom are the isospin doublet nucleon field 
^n(p), the meson fields which mediate the nuclear interactions 
(two isoscalar fields cr and co, two isovector fields p and n), 
and the photon field (A) that accounts for the electro-magnetic 
interaction. A general effective Lagrangian is constructed in 
terms of various currents or densities: (^TYxjj) where t e {1,t} 
(1 stands for the isoscalar quantities and the isospin Pauli ma¬ 
trix T corresponds to the isovector ones), and F denotes the 
Dirac matrices {l, 7 yi, 757 ^,(T^v}. Applying the standard vari¬ 
ations of the Lagrangian, one may obtain the field equations 
for nucleon, meson and photon fields (resp. the Dirac, Klein- 
Gordon and Proca equations), and the continuity equations 
leading to the Hamiltonian. With the creation and annihilation 
operators (c^,c) defined from stationary solutions of the Dirac 
equation, the effective Hamiltonian "K is formally expressed 
in the second quantized form as 


with the RHF model. It provides a unified and self-consistent 
description of both particle-hole {ph) and particle-particle 
ipp) channels, at the mean field level. 

B. Thermodynamics and statistical mechanics 

The thermodynamical properties of a system are calculated 
here in the canonical ensemble. For a statistical A-body sys¬ 
tem at finite temperature T, the equilibrium state is obtained 
from the variational principle applied to the grand canonical 
potential Qldl HUH, 

Q.{T,X}^F-A.N ^E-TS -AN, (3) 

where F is the free energy, S the entropy, E the total energy, 
and A the associated Lagrange multiplier. Namely, the varia¬ 
tion 

bQ = 0 (4) 

defines the density operator !D with trace equal to 1, and the 
grand partition function Z, respectively read as: 

D = Z“'exp{-^('H-dAO}, (5a) 

Z = Tr[exp{-/?C7T-dA())], (5b) 

where N is the particle number operator, and fi - T~^ . As 
it is conventional, the temperature T is given in energy units. 
For arbitrary operator O, the thermal average over the excited 
states populated at finite temperature is defined as: 

{O) = TxiDO), ( 6 ) 

where the trace is taken over all possible excited states. At 
finite temperature, the total energy, entropy and number of 
particles are therefore expressed as 

£ = (7T) = (7a) 

S = -<lnD) = -Tr(£)ln£)), (7b) 

A = (A() = Tr(2)A(). (7c) 

C. Extension to finite temperature 


ij ijkhip 

where Tij represents the Dirac kinetic energy, and the two- 
body terms Vijti correspond to different types of meson- (or 
photon-) nucleon couplings denoted by (p. In the mean-field 
approximation, the energy functional & is obtained by taking 
the expectation value of the Hamiltonian TY on a Slater de¬ 
terminant, where both the Hartree and Fock terms are consid¬ 
ered. The RHFB theory is deduced by incorporating 

the Bogoliubov transformation 


ru rv 


( 2 ) 


The FT-RHFB theory is a straightforward generalization 
of the RHFB theory that readily incorporates a statistical en¬ 
semble of excited states. At the finite-temperature mean field 
level, the density operator is approximated by d: 


where fa is the Fermi-Dirac distribution 


fa = {Na) = 


1 

1 + 


(9) 


Notice that, for T = 0 we have fa - 0 for all states a. The 
quasiparticle energy Ea is obtained as the solution of the FT- 
RHFB equations, see Eq. dnii. 
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In spherically symmetric systems the Dirac-Bogoliubov 
spinor can be written as 


^u,{r) ■■ 


1 


^v,{r) = - 


iGu^r) 

Fu,(r)o--r 

iGv^ir) ^ 
FvSr)a--r, 






(10a) 

(10b) 


where Gu and Fu correspond to the radial parts of the up¬ 
per and lower components, respectively, 'S^p^ir) is the spinor 

spherical harmonics, and Xi (t) is the isospinor. Here, the 
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subindex a —(a, m) j, m) contains the quantum numbers 
n (number of nodes of the upper component Gij(v)), I (orbital 
angular momentum), and j, m (total angular momentum and 
its projection to the z axis). 

The minimization of O with respect to D leads to the FT- 
RHFB equations. 




li{r,r') A(r,r') 
-A(r,r') h(r,r')j 


^v'(r') 


'A + E 0 ' 

'fuirf 

, 0 A-E, 

.fv(r), 


( 11 ) 


which are formally identical to the RHFB equations 

The Dirac Hamiltonian h{r,r') contains the kinetic energy h^, 

the direct local potential h^, and exchange non-local potential 


h^(r,r') = [a-p+(3M]6(r-r'), 

h^(r,r') = pr(r)r5 +So(r)+;825(r)]b(r-r'), 

„ (Yairy) YF(r,r')\ 
h {r,r ) = 

Ug(/-,/-') XF(r,r')} 


(12a) 

(12b) 

(12c) 


In the above expressions, the local self-energies 2s, 2o, and 
Sj- contain the contributions from the Hartree (direct) terms 
and the rearrangement terms, which depend directly on vari¬ 
ous local quasiparticle densities. They are the vector, scalar 
and tensor densities, respectively, 

a 

+[Gl^{r) + Fl^{r)]fa], (13a) 

a 

+[Gl^{r)-Fl^{r)]fa], (13b) 

a 

+2GuM)FuA>-)fa\- (13c) 

Notice that, at the limit T = 0, where - 0 Eqs. (I13ab - 
(I13cl) give back the usual expressions IHSll . The nonlocal 
self-energies Xq(f) and Yq(^f) come from the Fock (exchange) 
terms. 


+(8<i,Fuh)r^ab ("V’ ''')(8<l>Gui,)r'fb], (14a) 

- Z y^Vl[(8'pFyyA%^"(’”fYr'X8^FvM^ - fb) 

b 

+(8<pFUi)r^ab ("V’ Y r')(8^Fu,,)r'fb], (14b) 

Yg^yF) - Z ■^^«M[(^0Gv,).<«(m^;r,/)(g^Gv,).Kl -fb) 

b 

H84’Gub)rf^lhyFYr')(8,pGui,)r'fb\, (14c) 
Yf\Yr') = Z ^al?b[i8^Gvb)rA%l[im^-,ry)i8<^FvM^ -fb) 

b 

+ (8<l>Gui,)r^ah y<P’YF)(8<i>Fui,)r'fb]- (14d) 

In these expressions, denotes the isospin factors; dab for 
isoscalar channels and 2 - Sab for isovector channels, = 
2ji,+ l is the degeneracy number of the corresponding energy 
level, represents the coupling constants, and f^ab denote 
the multipole expansions of meson propagators. Eq. (fT4l i is a 
generalization at finite temperature of the expressions given in 
Ref. iS. 

Next, we consider the pairing field A(r, r') of Eq. (fTTI) . 

A„(r,/) = yj(r,r')^fc(r,/) (15) 

b 

with the pairing interaction V^'P and the pairing tensor k. If 
we take a finite-range pairing force, the pairing tensor k will 
read as 


Ka{r,r') = ]l{[GvAr)GuAF) + Fva(r)FuAr')] (16) 

+[GuAr)Gv,(r’) + Ec;„(r)Ev^,(/)]}(! - 2fa). 

Notice that the temperature dependence of the solution 
(Ea',fua’fVa) of the ET-RHEB Eq. (fTTI) comes implicitly 
through the quasiparticle densities pv,s,t, nonlocal potentials 
X{Y) and pairing tensor k. 

Eor practical evaluation of the pairing correlations, an av¬ 
erage pairing gap is introduced and defined as the ratio of the 
pairing energy over the pairing tensor. 


A = 


Tr(A/f) 
Tr/f ■ 


(17) 


This quantity is calculated for neutrons (A„) and protons (Ap) 
separately and it is discussed in detail in the next section. 

The total ET-RHEB energy E of the system is calcu¬ 
lated with the microscopic two-body center-of-mass correc¬ 
tion ii. The entropy S of the system can be evaluated from 


S(T) = [/«ln/„ + (l -/Jln(l -/„)], (18) 

a 


and the specific heat is defined by 


Cv(T) = 


^ dS(T) ^ 
dT U' 


(19) 


- Z ^aVl[<^gyVb)rf^abyFYr')i8^GvM^ - fb) 

b 


They correspond to the first and second derivative of the free 
energy F, respectively. 
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The integro-difFerential FT-RHFB Eq. (fTTl i is solved by us¬ 
ing a DWS basis with a radial cut off R = 26 fm. The 
numbers of positive and negative energy states in the basis ex¬ 
pansion for each s.p. angular momentum {I,]) are chosen to 
be 36 and 12, respectively. We have verified that this trun¬ 
cation scheme provides sufficient numerical accuracy for the 
description of weakly bound nuclei (Pb isotopes). 

III. RESULTS AND DISCUSSION 

In this section, we compare the predictions based on sev¬ 
eral models which are determined by the model Lagrangian 
as well as the pairing interaction. The latter interaction is ei¬ 
ther a finite-range Gogny DIS interaction lid]] or a density- 
dependent contact interaction (DDCI) of the form lld^ : 

V(r, /) = Vo 1 - Per) (1 - — ) - r'). (20) 

2 \ PO I 

Notice that the DDCI requires a regularisation scheme. We 
have considered in this study a simple cut-off scheme, defined 
to be 100 MeV in quasiparticle space, and the strength Vq is 
adjusted to reproduce the same average pairing gap as that 
obtained with the Gogny DIS interaction. These values of 
Vo will be given with the results hereafter. The DIS pairing 
interaction depends slightly on the mass: A general factor g 
is therefore introduced for its strength, as in Refs, if^ Isitl . 
We consider 6 different model Lagrangians which are given 
in Table H] Some of the bulk properties of nuclear matter de¬ 
termined by these Lagrangians are also shown. They are gen¬ 
erally compatible with the expected gross properties of finite 
nuclei, e.g., the binding energy E/A ^ 16 MeV, the satura¬ 
tion density po - 0.15 fm“^, the incompressibility modulus 
K e [220,280] MeV, and the symmetry energy J e [32,39] 
MeV. 

In the following discussion, the persistence and re-entrance 
of the pairing phenomenon will be commented and analysed. 
Let us briefly recall the unifying mechanism which is at pla y 
in these various phenomena, recently discussed in Ref. lU^ . 
Due to thermal excitations, s.p. states above the Fermi en¬ 
ergy can be slightly populated while states below the Fermi 
energy can be partially depleted. This occurs if the involved 
new states are not too far in energy from the last occupied 
state, but it should also be not too close, otherwise, these states 
would already participate to pairing correlations at zero tem¬ 
perature. The typical shell gap should be around 2 MeV. The 
participation of these states at finite temperature gives rise ei¬ 
ther to the persistence of pairing correlations slightly above 
the usual critical temperature for nuclei which are already su¬ 
perfluid at zero temperature, or to pairing re-entrance at finite 
temperature for nuclei which have weak or no pairing at zero 
temperature. The best nuclei, in which such a phenomenon 
is expected, are those close to the drip line, as well as those 
located at a subshell closure as shown in this section. 

In the following subsections, we evaluate the influence of 
the model on the pairing properties in hot finite nuclei, tak¬ 
ing as an example for testing pairing correlations. The 
ratio of the critical temperature over the average pairing gap 



T (MeV) T (MeV) 

FIG. 1. (Color online) The neutron pairing gaps in '^“^Sn as a func¬ 
tion of temperature, calculated with PK02 (black) and DD-ME2 
(blue) corresponding to different In the pairing channel, we 

compare finite-range DIS (filled circles) and contact DDCI (empty 
circles) forces. The analytical results (dashed lines) are also shown. 
Notice that here, we took the factor g = 1 for the Gogny DIS force, 
and for the DDCI pairing, we have taken the following values for Vq 
(in MeV.fm^^): 335 (PK02), and 342 (DD-ME2) for the RH(F)B 
and BCS calculations. Eor the DDCIx force used in BCS, we took 
for Vq (in MeV.fm-^): 526 (PK02) and 539 (DD-ME2). 

at zero temperature Tc/A„{Q) is also studied as a function of 
the model, varying either the Lagrangian or the type of pair¬ 
ing interaction. Linally, a more systematic study of the critical 
temperature is performed on a set of semi-magic nuclei. 

A. Study of the ratio Tc/AnlO) 

Since pairing correlations are active only around the Lermi 
level, the ratio Tc/A„(Q) is expected to be modified by the ef¬ 
fective mass which influences the s.p. level spacing to a large 
extent. Notice that here, the effective mass corresponds to the 
non-relativistic one instead of the mantity — M + Y^s 
that is named as the Dirac mass 1^ . It can indeed be 
shown that, in the weak coupling limit of the BCS approxima¬ 
tion, the average pairing gap at the Fermi surface Af can be 
expressed as 

Af «;2ef exp[2/(VfVpair)], (21) 

where ep is the Fermi energy, Np - nfpkp/jp' is the average 
density of state in uniform matter at the Fermi energy, and kp 
denotes the Fermi momentum, and Vpaij is a constant pairing 
interaction. It is clear from Eq. (i2Ti) that the pairing gap Ap 
is quite sensitive to the effective mass at the Lermi energy ep. 
Eq. (i2Tl i is obtained in infinite matter and it provides only a 
qualitative understanding of the relation between the pairing 
force strength and the effective mass. In the following, we 
present a quantitatively precise analysis of the correlation be¬ 
tween the critical temperature and the non-relativistic effective 
mass in finite nuclei. 

Figureina) displays the evolution of the neutron pairing gap 
as a function of temperature for *^^Sn, a good candidate for 
studying pairing correlations. We compare two Lagrangians, 
PK02 and DD-ME2 (see Table|I]i both with two kinds of pair¬ 
ing interactions: the finite-range Gogny DIS force IdUl and 
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TABLE I. Bulk properties of symmetric nuclear matter at the saturation point: density po (fm ^). binding energy Eg!A (MeV), compression 
modulus K (MeV), symmetry energy J (MeV) and non-relativistic effective masses (M) predicted by selected RHF and RH models. The 
non-relativistic effective masses in neutron matter are also listed. 


Model 

Interaction 

Ref. 


symmetric matter 


neutron matter 

Po 

Eg/A 

K 

J 

^NR 

MnrM 

AInr^^) 


PKAl 


0.160 

-15.83 

229.96 

36.02 

0.68 

0.68 

0.70 

RHF 

PKOl 


0.152 

-16.00 

250.28 

34.37 

0.75 

0.73 

0.76 


PK02 


0.151 

-16.03 

249.53 

32.49 

0.76 

0.75 

0.77 


DD-ME2 


0.152 

-16.14 

250.97 

32.31 

0.65 

0.64 

0.70 

RH 

PKlr 


0.148 

-16.27 

283.68 

37.83 

0.68 

0.64 

0.72 


NL3* 


0.150 

-16.30 

258.56 

38.70 

0.67 

0.63 

0.72 


the contact force DDCI lld^ . The two Lagrangians PK02 and 
DD-ME2 mostly differ by their non-relativistic effective mass 
(see Table 0, and it is observed that the average pair¬ 
ing gap at zero temperature A„ scales with as expected 
from the weak coupling expression (l2Tl i. Comparing the dif¬ 
ferent types of pairing interaction (finite- or zero-range) for 
the same Lagrangian, it is observed in Fig.[T]that the vanishing 
of pairing correlations at finite temperature slightly depends 
on the type of pairing force, namely, with the zero-range pair¬ 
ing interaction, the critical temperature Tc is slightly lower 
than with the finite-range interaction. From Fig. [1] the ratio 
TrlA„(0) is obtained as 0.60 for DIS, and 0.57 for DDCI. In 
Ref. 11^ . it was shown that the dependence of the pairing gap 
on the state around the Fermi energy is qualitatively different 
for contact and finite-range interactions. In addition, with a 
finite-range interaction, the pairing gap and the pairing ten¬ 
sor have non-local components which cannot be simply ab¬ 
sorbed in the DDCI. The slight increase of the critical tem¬ 
perature with the DIS interaction is therefore an effect of the 
finite-range nature of the interaction, dispersing the pairing 
effects among more s.p. states. For comparison, and taking 
Tc = 0.60A„(0) (finite-range pairing force) and Tc - 0.57A„(0) 
(zero-range pairing forces), the analytical relation if^ . 


A„(7’) = A„(0) 


1 - 



( 22 ) 


readjusted Vq in the RHF-BCS calculation at zero tempera¬ 
ture to obtain the same pairing gap as the RHFB prediction, 
leading to the DDCIx interaction in Fig.[TJb). Applying such a 
simple modification of the parameter Vq in the tempera¬ 
ture dependence of the pairing gap predicted by the FT-RHFB 
and FT-RHF-BCS frameworks are almost undistinguishable. 
However, this is not always true and as we will see that the 
above simple renormalization of the pairing strength will not 
work towards the drip line where the coupling to continuum 
states becomes more and more important. 

The critical temperatures is also expected to depend on 
the non-relativistic effective mass, see for instance Refs, ll^ 
The dependence of the ratio Tc I A„{0) on the effec¬ 
tive mass is, however, not very well known. In Fig. |2] we 
plot the ratio Tc I A„(0) as a function of in '^‘^Sn, and 
we compare the predictions of two pairing interactions (fi¬ 
nite versus zero range). The effective mass is obtained 
with 14 CDFs (the 6 CDFs given in Table U completed with 
8 other CDFs: PK03 JH, PKDD lH, DD-MFl IHl], DD- 
MFb JH, TW99 JH, PKl iH, NL3 H and TMl iH). 
It is found that the ratio Tc I A„{Q) does not depend much on 
M*^^. Similar analyses were also carried out to check the rela¬ 
tions with the incompressibility modulus K and the symmetry 
energy J, but no evidence of correlation is found. The small 
fluctuations observed in Fig. |2] are therefore mostly shell ef¬ 
fects. In conclusion, the critical temperature in stable nuclei 


where m — 3.32, is plotted for the DD-MF2 model [dashed 
lines in Fig.[TJa)]. There is almost no difference between the 
analytical model, i.e., Fq. (l22li . and the numerical calculations 
for stable nuclei like ^^^Sn. 

Figure [Ttb) shows the neutron pairing gap calculated with 
the RHF theory plus BCS pairing at finite temperature (FT- 
RHF-BCS), using the same DDCI interaction as in the FT- 
RHFB calculations shown in Fig. [Ha) and modified one (DD¬ 
CIx) with enhanced pairing strength parameter Vq (see cap¬ 
tion for more details). It is found that, with the same DDCI 
pairing interaction, the neutron pairing gaps determined by 
the BCS method are reduced to about half of Bogoliubov re¬ 
sults. Such distinct difference is due to the fact that the off- 
diagonal couplings, which account for about half of the pair¬ 
ing correlations, are absent in the BCS pairing. As a result, the 
strength parameter Vq in BCS calculations is usually larger 
than in HFB, see for instance Ref. if^ . We have therefore 



FIG. 2. (Color online) The ratios TdA„(0) in *^^Sn as a function 
of the non-relativistic effective mass calculated using the FT- 

RH(F)B theory with 14 parameter sets. In the pairing channel the 
finite-range DIS and the contact force DDCI are employed. 
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FIG. 3. (Color online) The critical temperature Tc (left axes, blue 
circles) and the occupation number of continuum states Neon (right 
axes, black squares) at zero temperature in *^^Sn as a function of 
neutron pairing gap calculated from PKOl FT-RHFB with Gogny 
DIS (a) and DDCI (b) pairing forces. The dashed lines correspond 
to the relation Tc = 0.60(0.57)A„(0). The strength parameter of the 
pairing interaction is varied to obtain different values of the pairing 
gap at zero temperature. 


scales very well with the average pairing gap at zero tempera¬ 
ture, the shell effects contributing to a dispersion of less than 
4%. 

In order to check whether the ratio 7’c/A„(0) is still con¬ 
stant for larger pairing strength, we have artificially varied 
the pairing interaction strength and correlated the critical tem¬ 
perature Tc with the neutron pairing gap A„(0) at zero tem¬ 
perature. Figures |3ja) and [3b) respectively show the re¬ 
sults calculated with the pairing interactions DIS and DDCI 
using the PKOl Lagrangian, where the average pairing gap 
A„(0) goes from 100 keV up to about 3.5 MeV. At low pair¬ 
ing gap, the ratios Tc I A„{Q) are consistent with the analytic 
ones (dashed lines) as expected, until small deviations emerge 
beyond A„(0) ~ 2.5 MeV, where the continuum contributions 
become sizable (filled squares). These results are consistent 
with previous findings l!^ based on a separable version of 
the Gogny DIS pairing force ll^ and RH Lagrangian with 
PC-PKl point coupling l!^ . A simple calculation in infinite 
matter with a contact pairing interaction gives a correction to 
the ratio ll^ 


^»0.57|i-^A(0)T (23) 

where a>D represents the pairing window. According to 
Eq. (i23 l. the next-to-leading order correction has a negative 
sign, at variance with our results in '^‘^Sn, see Fig. [3] Notice 
that for the DDCI pairing interaction, the pairing window is 
about 100 MeV. The correction to the linear approximation of 
Eq. ( l23 l is therefore very small; for the maximal pairing gap 
considered in this work (A(0) « 4.0 MeV), the correction rep¬ 
resents no more than a few percent of the linear leading term. 

In Fig. 13 the increase of the critical temperature for A„(0) > 
2.0 MeV reveals an enhancement of the thermal pairing cor¬ 
relations, as well as an the important role played by the con¬ 
tinuum states. We remind that the next-to-the-leading order 
correction appearing in the simple expression in uniform mat¬ 



FIG. 4. Contributions, in *^^Sn, of the continuum states to the pair¬ 
ing number Neon (a) and to the neutron number Neon (b), normalized 
to their value at zero temperature, as a function of the temperature 
T/O.bAfO). The results correspond to PKOl Lagrangian and DIS 
pairing interaction. 


ter (|23 is negative. However Fig. [3] shows a continuous en¬ 
hanced continuum effect with respect to the pairing gap A„(0). 
To better stress the increasing contribution of the continuum 
states and their role in the pairing correlations, we have intro¬ 
duced the following two quantities: the cumulative occupation 
number of the neutron continuum states Neon, 

Neon — 2 r Anr^pa{r)dr, (24) 

a,ea>Q 

and the cumulative pairing-occupation number of the neutron 
continuum states Neon, 

Neon ^ Xj f I^aiddr. (25) 

a,€c>0'^ 

Here the continuum states a are determined in the T = 0 
canonical basis, for simplicity Elil, with s.p. energy ea 
above the continuum threshold. The increasing role of the 
continuum states at finite temperature is illustrated in Fig. |4| 
where the contributions to the pairing number Neon [plot (a)] 
and to the neutron number Neon [plot (b)] from the continuum 
states in ^^^Sn, normalized to their values at zero temperature, 
are shown. Results with a weak pairing (A„(0) =1.3 MeV) and 
with a stronger pairing (A„(0) = 3.0 MeV) are compared. Even 
if in the latter case the pairing is slightly larger than the ex¬ 
pected value in finite nuclei, its inclusion in our analysis helps 
to understand the role of the continuum states. For the weak 
pairing case, the continuum effects are very small, see Fig. |3 
and both Neon and Neon drop to zero at the expected value 
Te - 0.6A(0), see Fig. |4| For the strong pairing case, a clear 
correlation is observed in Fig. |4| between the increase of the 
occupation of the continuum states at finite temperature Neon 
[see plot (b)] and the persistence of the pairing numbers Neon 
[see plot (a)]. The persistence of pairing correlations below 
the critical temperature modifies also the critical temperature. 
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FIG. 5. (Color online) Comparison of A(0) (black circles), 0.60A(0) 
(green curves) and Tc (blue circles) in the even-even Ni, Sn, Pb iso¬ 
topes (left panels) and N = 50, 82, 126 isotones (right panels), cal¬ 
culated in FT-RHFB with PKAl and the Gogny pairing interaction 
DIS. 


and it is observed in Fig. Ufa) a larger value of the quantity 
T /0.6A(0) where pairing correlations in the continuum space 
drops to zero in the case of strong pairing compared to the 
weak one. Coming back to Fig. [3] we now understand bet¬ 
ter the correlation between the slight deviations of Tc from 
the analytical behaviors and enhanced continuum effects for 
A(0) > 2.5 MeV. Since the presence of resonance states in the 
continuum is a typical feature of finite systems, the increase 
observed in Fig. [3] which is different from the prediction of 
Eq. (EDI, is then expected only in finite systems. Anticipating 
the results shown in the next figures, a similar enhancement 
of pairing correlations in other nuclei will also be observed, 
revealing here also the role of the resonant states. 


B. Evolution of the critical temperature in isotopic and 
isotonic chains 

In this subsection, we perform a systematic analysis of the 
evolution of the pairing gaps and critical temperature along 
isotopic and isotonic chains of semi-magic nuclei. Through 
this extensive analysis, we have access to various s.p. config¬ 
urations going from stable nuclei towards weakly bound drip 
line nuclei, and we probe the pairing correlations inside vari¬ 
ous major shells. We consider three models, PKAl and PKOl 
(RHF) and DD-ME2 (RH) which have different symmetry en¬ 
ergies and non-relativistic effective masses, see TableH] Antic¬ 
ipating our results, we will see that these models lead to rather 
different predictions for the pairing gap A. 

Our results are shown for PKAl (Pig. |3, PKOl (Pig. [b]) 
and DD-ME2 (Pig. [T]) models, with the pairing channel de¬ 
scribed by the DIS interaction. In Pigs.|5]l2l we have repre¬ 
sented isotopic (Ni, Sn and Pb) and isotonic {N - 50, 82, and 
126) average pairing gap evolution as a function of N (for iso¬ 
topes) and Z (for isotones). The isotopic and isotonic chains 
are bounded by the drip lines predicted by each of the con¬ 
sidered models and determined by the two-nucleon separation 
energy. These drip lines are consistent, within a few units un¬ 
certainties, with predictions given by other models obtained 
with Skyrme forces Gogny forces IiIItI and RH 

Pagrangians if^ Ft^ . We have calculated the average pair¬ 
ing gap at zero temperature (filled black circles) defined from 
Eq. (O, and compared the calculated critical temperature 
(filled blue circles) with the approximate relation, i.e., 0.6A(0) 
(green curves). The arch structure of the results shown in 
Pigs. |5]l2] reflect the presence of magic numbers where pair¬ 
ing correlations completely vanish. 

The PKAl model (Pig.|5j is the most complete RHP version 
of the CDP theory. It contains the p-N Porentz tensor coup ling 
which is known to enhance the spin-orbit splitting 04311^1^ : 
in many cases the subshell structure is found to be closer to 
the experimental data than those predicted by other models 
without the p-N Porentz tensor coupling, such as the RH ap¬ 
proaches shown in Pig.|7] These subshell structures are clearly 
visible in Pigs.|5]l7]since they induce a partial quenching of the 
pairing gap for the associated submagic numbers. Going to¬ 
wards the drip lines, a reduction of the pairing gaps is often 
observed, revealing the presence of closed-shell nuclei at or 
near the drip lines. Por the neutron drip line, it is the case of 
Sn and Pb isotopes, and for the proton drip line, it is observed 
for N = 50. 

We first discuss the pairing properties of finite nuclei at 
zero temperature, which are influenced by the underlying s.p. 
structure around the Permi energy. Por the Ni isotopes, a sub¬ 
shell closure at A = 40 is predicted with PKAl and PKOl Pa¬ 
grangians, as expected from experiments iIt^ . while DD-ME2 
shows a more pronounced shell closure. Por neutron rich Ni 
isotopes, PKAl indicates another subshell closure at A = 62 
which is not seen with PKOl or DD-ME2. It is, however, 
beyond the present experimental limits. Por the Sn isotopes, 
a decrease of the pairing gap induced by subshell closure is 
observed at A = 64 with PKOl and DD-ME2 Pagrangians, 
but not with PKAl. Por Pb, a small decrease of the pairing 























8 



FIG. 6. (Color online) The same as in Fig. [5] but calculated in the 
FT-RHFB theory with the effective interaction PKOl. 


gap is observed nt N — 138 with PKAl and nt N — 146 with 
PKOl and DD-ME2. On the isotonic side, we observe a sub¬ 
shell closure at Z = 40 for N = 50 with PKAl Lagrangian, 
but not with PKOl or DD-ME2. Eor N - 82 isotones, PKAl 
predicts a well marked shell closure at Z = 40 and a subshell 
closure at Z = 34 and Z = 64, while at Z = 58 PKOl predicts 
a shell closure and DD-ME2 only a subshell closure. Einally, 
for N - 126, PKOl and DD-ME2 indicate a reduction of the 
pairing gap at Z = 92, which is not confirmed by experimen¬ 
tal data and is not present with PKAl Lagrangian. We 
will see below that these structures can have an impact on the 
thermal properties. 

Turning now to the thermal properties of these isotopes and 
isotones, the comparison of the calculated critical temperature 
Tc and the approximate relation 0.60A(0) shown in Eigs. |5} 
|7] exhibits some interesting features. The critical tempera¬ 
ture Tc and the approximate relation 0.60A(0) are identical 
in most cases with some exceptions. In heavy nuclei (Pb and 
N - 126), there are no strong deviations between these two 
quantities, but they are however more marked in lighter nu¬ 
clei. Moreover, the cases where the exact and the approximate 



FIG. 7. (Color online) The same as in Fig. [5] but calculated in the 
FT-RHB theory with the effective interaction DD-ME2. 


values of Tc differ are correlated with either the presence of a 
subshell closure, or with the proximity of the drip lines. In 
the case of subshell closure, the effect of the temperature is to 
’’wash out” the decrease of the pairing correlations. This can 
be understood as the consequence of the thermal occupation 
probabilities which overcome small shell gaps. Close to the 
neutron drip lines, the more pronounced effects are observed 
in Ni and Sn isotopes. This is also due to the thermal occupa¬ 
tion of close-by resonant states as discussed in Refs. 0141^ . 
In the non-relativistic Skyrme Hartree-Eock plus BCS (SHE- 
BCS) approach, an enhancement of the critical temperature 
was found in using SkT6 1^ . We do not conhrm 

this enhancement in '“^^Sn with the models used in this work. 
However, it is interesting to notice that the origin of such an 
enhancement found in Ref. is also related to the existence 
of a subshell closure. 

Let us finish this subsection with some general remarks 
concerning the nuclei which do not manifest any enhancement 
of the critical temperature. Eor the Pb isotopes, as shown in 
Pigs. 137] we have not observed any marked enhancement of 
the critical temperature near the drip line as in the case of 
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T (MeV) T (MeV) 

FIG. 8. (Color online) The neutron pairing gaps in ®^Ni as a function 
of temperature, calculated by FT-RHFB with Gogny DIS and DDCI 
pairing forces, and by FT-RHF-BCS with DDCI pairing force. To 
evaluate the persistence provoked by the subshell, the analytical re¬ 
sults are also shown. The pairing strength Vq (in MeV.fm”^) is fixed 
to be 326 (DDCI) and 537 (DDCIx). 


Sn isotopes. Comparing Pb and Sn, since the pairing gap for 
these isotopes is decreasing near the drip line, one could have 
expected to observe an enhancement of the critical tempera¬ 
ture in Pb as it is observed in Sn. For instance, the last oc¬ 
cupied state in the drip line nucleus 266pb is indeed found to 
be well bound ^ “^.0 MeV), and the lowest s.p. reso¬ 

nance 62 / 111/2 found to be above 1.5 MeV. There is therefore 
a rather strong gap in the neutron rich Pb isotopes {N = 184) 
which prevents the large coupling to the continuum. For the 
isotonic chains, we do not find any pairing persistence phe¬ 
nomenon around the drip line. This can be related to the well 
developed shell closures at Z = 50,82 and 92 for proton-rich 
nuclei, which quench the coupling to continuum states, as it 
would have been expected in such exotic nuclei. In addition, 
the coupling to the continuum is weaker for protons since the 
Coulomb barrier tends to localize the proton density in the 
nuclear interior if^ . For these reasons, the persistence phe¬ 
nomenon is strongly quenched in the proton channel. 


C. Pairing persistence in ®^Ni and *^^Sn 

In this subsection, we analyze in more details the tempera¬ 
ture dependence of the pairing gap for two representative nu¬ 
clei, ®^Ni and *^^Sn. On the one hand, ®^Ni is an isotope 
which is slightly more neutron rich than the five stable iso¬ 
topes As shown in the prev ious subsection, ®^Ni is 

considered as a subshell isotope 11 ^ 1 ^ 1 ^ . and as a conse¬ 
quence, the pairing gap at zero temperature is either reduced 
or strongly quenched depending on the model, see Figs. EE 
On the other hand, *^^Sn is a very neutron-rich isotope at or 
close to the neutron drip line, where the continuum effects are 
expected to be remarkable dlllllllll. However, since 
*^^Sn is close to the potentially doubly magic *^®Sn (Z=50, 
A=126), a gap is expected to be present in the s.p. struc¬ 
ture between bound and resonant states. These two nuclei are 
therefore representative of quantum systems for which pair¬ 
ing at zero temperature is weakened by the presence of a gap 


above the Fermi energy and where a finite amount of tempera¬ 
ture allows to overcome the gap and provokes an enhancement 
of pairing correlations, giving rise to pairing persistence. 

We first show the temperature dependence of the pairing 
gap for ®^Ni in Fig. 0a), calculated with the FT-RHFB model 
where we consider the PKOl Lagrangian in the mean field 
channel and either Gogny DIS or DDCI interaction in the 
pairing channel. The analytical model is also shown for ref¬ 
erence. It is found that the predictions for A„{T) do not prac¬ 
tically depend on the pairing force. The critical temperature 
predicted by the FT-RHFB approach is increased with respect 
to the reference analytical model: the FT-RHFB pairing gaps 
vanish around T - 0.90 MeV, which is 0.25 MeV higher than 
the expected value (0.65 MeV). The pairing gap predicted by 
FT-RHF-BCS is shown in Fig. 0b). Surprisingly, the pair¬ 
ing gap is zero if the same DDCI pairing interaction is used. 
An increase of the pairing strength Vq is therefore necessary. 
It is also interesting to observe that the DDCIx pairing inter¬ 
action, where Vo is increased to match with the zero temper¬ 
ature pairing gap obtained with FT-RHFB case, reproduces 
almost exactly the temperature dependence of the FT-RHFB 
case and predicts as well an increase of the critical tempera¬ 
ture with respect to the analytical model. The nucleus ®^Ni is 
a typical example of pairing persistence at finite temperature 
in closed subshell {N = 40) nuclei. We hereafter name this 
phenomenon type I pairing persistence. Other examples of 
similar behavior are: ‘^‘’Ni (PKAl), (PKOl, DD-ME2), 
220 pb (PKAl), ^^‘’Pb (DD-ME2) for the neutron pairing gap, 
and ^OZr (PKAl, PKOl and DD-ME2), '^Oce (DD-ME2), 
146,190 qj (PKAl) for the proton pairing gap. 

We turn now to the analysis of the results in *^^Sn. As 
stated above, this is a nucleus where pairing correlations are 
slightly weakened due to the proximity of shell closure. A 
small amount of temperature is expected to reorganize the 
level occupancy around the Eermi energy, opening more space 
below the Eermi energy, and producing a non-zero occupancy 
of the first levels above the Eermi energy which are in the 
continuum. Most of the occupied states in the continuum are 



T (MeV) T (MeV) 

FIG. 9. (Color online) The neutron pairing gaps in *^^Sn as a func¬ 
tion of temperature, calculated with PKOl and NL3*, using Gogny 
DIS and DDCI pairing forces. The results of the analytical model are 
also shown. The pairing strength Vq (in MeV.fm^^) is: 333 (DDCI 
with PKOl), 317 (DDCI with NL3*) and 596 (DDCIx with PKOl 
and BCS framework). 
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resonance states, but it is interesting to notice that a small 
number of them are also non-resonant states JUt]. Without 
the participation of these non-resonant states, the asymptotic 
behavior of the density would be ill-defined and present an 
unexpected gas component, as it is also observed in the BCS 
theory 0. The Bogoliubov transformation couples all states 
in a sub-(/, 7 ) space and a truncation among these states breaks 
the unitarity of this transformation. To avoid the presence of 
non-physical gas component in the density profile, it is there¬ 
fore important to couple all states in the continuum within the 
Bogoliubov transformation. 

Figure |9ja) displays the evolution of the neutron pairing 
gap as a function of temperature in ^^^Sn with different La- 
grangians and pairing forces. The results are very weakly 
affected by the choice of the pairing interaction. The effect 
induced by the choice of Lagrangian is also very small. We 
find a systematic increase of the critical temperature {Tc ~ 
0.65 - 0.70 MeV) with respect to that expected from the an¬ 
alytical relation {Tc ~ 0.47 MeV), independently of the con¬ 
sidered model. In addition, for temperature T > 0.2 MeV, an 
increase of the pairing gap is observed, which can also be re¬ 
lated to the thermally induced contribution of the continuum 
states. We compare these results to the ones obtained with the 
FT-RHF-BCS framework shown in Fig.|9tb). As already dis¬ 
cussed, the DDCI interaction predicts a reduced pairing gap 
in BCS compared to RHFB. On the other side, the DDCIx 
interaction where the pairing strength is increased to match 
the T - 0 predictions of RHFB leads to an overestimation of 
the critical temperature compared to the FT-RHFB case. It 
shows that, in this case of dripline nucleus, the RHFB calcula¬ 
tion cannot be simply reproduced by a BCS calculation where 
the pairing strength is just increased. Since the coupling to 
the continuum plays a dominant role in the persistence phe¬ 
nomenon in *^^Sn, we hereafter call it type II phenomenon. 
From our results, it is also expected to occur in Ni and Sn 
neutron-rich nuclei, Ni (V > 54 ~ 60) and Sn {N > 100). 


D. Entropy and specific heat 

We now focus on the entropy S and specific heat Cy which 
are the first and second derivatives of the free energy F 
with respect to the temperature, and thus sensitive to thermal 
changes of the ground state, see for instance Ref. ll^ and 
Refs, therein. To test the sensitivity of these quantities to the 
choice of different models, we select two Sn isotopes, ^^^'Sn 
and ^^'^Sn. The former is a good example of a stable well- 
paired nucleus, while the latter is yet far from actual nuclear 
experiments but represents an extreme case with large isospin 
asymmetry. 

In Fig. [To] are shown the entropy S and the specific heat Cy 
as functions of the temperature calculated by the RHF func¬ 
tionals PKOl and PK02, the RH ones with the non-linear 
self-couplings PKlr and NL3*, and the RH one with density- 
dependent meson-nucleon couplings DD-ME2. In the pairing 
channel, the value of the scaling factory (see first paragraph of 
section EDl, is slightly modified to give identical pairing gaps 
at zero temperature for the different models. In Fig. fTOl' a). the 
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FIG. 10. (Color online) Entropy and specific heat in ^^^Sn as a func¬ 
tion of temperature, calculated using the FT-RH(F)B and FT-RH(F) 
theories. 


entropy calculated with the pairing correlations [labelled by 
FT-RH(F)B] is compared to that neglecting the pairing cor¬ 
relations [labelled by FT-RH(F)]. At low temperature, if the 
pairing effects are ignored, the entropy is found to be largely 
model-dependent, i.e., the model with smaller non-relativistic 
effective mass (see Table |I]l, which leads to larger s.p. level 
spacing on the average, presents smaller entropy. As the tem¬ 
perature increases, and also as the pairing correlations are 
switched on, the entropy becomes less model-dependent. 

In fact, at low temperature or without pairing correlations, 
the entropy is largely determined by the few states around 
the Fermi energy, and the number of the involved states is 
essentially determined by the detailed s.p. spectrum which 
depends on the models, therefore leading to model-dependent 
entropy. Both temperature and pairing correlations can dis¬ 
perse the particle over the states beyond the Fermi level. As 
the temperature increases, and/or as the pairing correlations 
are enhanced, more s.p. states will get involved to contribute 
to the entropy, and the average properties such as the den¬ 
sity of states will become dominant, instead of a few states 
in the FT-RH(F) cases. Compared to distinctly different s.p. 
spectra around the Fermi surface, the dispersions of the non- 
relativistic effective masses (see Table HI between the models 
are less remarkable. Even though, in Eig. [TOl'a) it is clearly 
shown that the ET-RH(E)B results are grouped by the values 
of the effective masses when T > 1 MeV which correspond 
to different average densities of states. As expected, the ef¬ 
fect of the pairing correlations is clearly visible below the 
critical temperature {Tc ~ 0.8 MeV), inducing a strong reduc¬ 
tion of the entropy [see Eig. [TOl' all and singular behaviors of 
the specific heat around the critical temperature as shown in 
Eig. lTOl bl. Just above the critical temperature, we can notice 
that the specific heat is not linear in T, as expected from the 
Eermi gas model and the linear dependence seems 

to be found at slightly larger temperature (T > 1.5 MeV). The 
non-linearity of the specific heat around Tc might be related 
to shell effects. 

The results thus clearly show that the pairing correlations 
contribute to the s-shaped behavior of the specific heat, as 
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FIG. 11. (Color online) The same as in Fig.[T0l but for *®°Sn. 

it has already been inferred from the anal ysis of thermal ex¬ 
cited nuclei in laboratory experiments 11 a lift 1^ . A realis¬ 
tic description of the smooth s-shaped behavior in finite nu¬ 
clei requires a more elaborated modellingjincluding for in¬ 
stance particle number projection O, ISTI - I^ . It is however 
shown in Ref. ll^ that the smooth i-shaped behavior may 
be even washed out in some rare earth nuclei. The results 
presented in Fig. [10] should not be compared directly to the 
semi-experimental data. 

The situation for the neutron-rich nucleus ^®**Sn is more 
complex as shown in Fig. [TT] The model dependence of en¬ 
tropy at low temperature is reduced compared to the case of 
Only PK02 predictions differ from the other models. 
The predictions for the critical temperature, associated to the 
discontinuity in the specific heat of Fig. [TT] vary among the 
modellings to a much larger extent than what was found in 
I20sn jg jocated in the region where pairing persistence 

is expected to appear, see Figs. EH Since this phenomenon 
is strongly related to the position of resonance states in the 
continuum, we expect to observe deviations among models 
predicting different positions of these states. This model de¬ 
pendence therefore reveals our lack of knowledge in extrap¬ 
olating Lagrangians which have been adjusted for less exotic 
nuclei. 


IV. SUMMARY AND CONCLUSION 


range nature of the pairing force, while generating different 
state-dependence pairing gaps, causes only small differences 
in our results. We have described the pairing persistence in 
two kinds of situations: nuclei at subshell closure (type I), and 
nuclei strongly coupled to continuum states which are close to 
the drip line (type II). We have observed that, while a refitting 
of the pairing strength could match the FT-RHF-BCS with the 
FT-RHFB predictions for the pairing gap in the case I, it is 
no longer the true in the case II. This is due to the participa¬ 
tion of the continuum states in the second case which involve 
coupling of different nature. 

We have also analyzed the influence of the interaction on 
the thermal response. The results show clearly that the pair¬ 
ing correlations contribute to the i-shaped behavior of the spe¬ 
cific heat curve, and help to wash out the model dependence. 
For stable nuclei the model deviations, to some extent, can 
be traced back to the effective mass, since the level structure 
only weakly depends on the choice of the CDF. The situation 
for exotic nuclei is more complex since it is related to our lack 
of knowledge in very exotic nuclei, and the pairing persistence 
would have large effects on their thermal property. 

In conclusion, we have illustrated the richness and com¬ 
plexity of pairing correlations at finite temperature and in fi¬ 
nite systems within the first FT-RHFB calculation. The dis¬ 
cussion of correlations beyond mean-field, induced for in¬ 
stance by particle number projection, is not addressed in this 
work. It is however expected that the particle number pro¬ 
jection will contribute to increase the pairing correlations in 
the case where they are weak ll^ . like in the pairing persis¬ 
tence phenomenon discussed in this paper. In future work, 
a more quantitative calculation will be necessary to estimate 
how strong are these additional correlations and how they 
modify the results presented in this work. Another interesting 
perspective which is suggested by this work is the possibil¬ 
ity that similar phenomena can be observed in other domains 
of physics. For instance, it was studied whether cold atoms 
in a double potential could demonstrate pairing persistence as 
well ll^ . Finally, the application of this formalism for the 
prediction of temperature evolution of pairing properties in 
the crust of neutron stars ITTl 193] will be performed in the 
near future. There, the thermal modification of pairing cor¬ 
relations could have a large impact on the thermal relaxation 
of the crust iH, and could be observed during the quiescent 
period of low mass X-ray transients. 


In this work, we have developed the first FT-RHFB the¬ 
ory for spherical nuclei. The self-consistent FT-RHFB equa¬ 
tions are solved by using a DWS basis which provides an ap¬ 
propriate asymptotic behavior for the continuum states. We 
have performed systematic FT-RHFB calculations for both 
stable and weakly bound nuclei and discussed their thermal 
properties. The influence of the pairing interaction on the 
pairing phase transition is evaluated. It is found that the 
critical temperature for a pairing transition generally follows 
the rule Tc — 0.60A(0) with a finite-range pairing force and 
Tc — 0.57 A(0) with a contact pairing force. The finite- or zero- 
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